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Abstract. Model of biological effects of ionizing particles, especially of protons and 
other ions, is proposed. The model is based on distinguishing the single-particle and 
collective effects of the underlying radiobiological mechanism. The probabilities of 
individual particles to form severe damages to DNA, their synergetic or saturation 
combinations, and the effect of cellular repair system are taken into account. The 
model enables to describe linear, parabolic and more complex curves, including 
those exhibiting low-dose hypersensitivity phenomena, in a systematic way. Global 
shape as well as detailed structure of survival curves might be represented, which is 
crucial if different fractionation schemes in radiotherapy should be assessed precisely. 
Experimental cell-survival data for inactivation of V79 cells by low-energy protons have 
been analyzed and corresponding detailed characteristics of the inactivation mechanism 
have been derived for this case. 
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1. Introduction 



Models estimating the fraction of surviving cells are of utmost importance for 
applications of ionizing radiation in radiotherapy, because both tumor control and 
normal tissue complication probabilities (TCP, NTCP) are given primarily by cell 
survival in the corresponding regions. In the case of hadron radiotherapy, based on 
irradiating by accelerated protons and light ions, the biological effects are not given 
by applied dose only, but depend significantly on other physical characteristics of these 
particles. Their relative biological effectiveness ( RBE) has been shown to vary hugely 



with ion kind and energy or LET value; see, e.g., ((Belli et a,l Ill998l 1200(1 iFolkard et al 



1996 



Wevrather et a/l ll999V 



Several phenomenological model approaches aimed at estimating the fraction of 
surviving cells have been proposed so far for treatment planning procedur es, based e.g 



on a similarity in LET characteristics between carbon and neutron beams (|Kanai et al 



1999) or using the linear-quadr atic (LQ) forma lism with its a and (3 coefficients 
depending on LET and ion kind ((Belli et al 1119971 ) . A systematic m odel approach has 
been represented by the local effect model fLEM) of Scholz and Kraft ((Scholz and Kraft 
1994L Il995l IScholz et al II1997L Kraft et anll99 9). in which the differences in biological 
effects of diverse ions have been related to their track structures. However, in any of 
these approaches the underlying physically-chemical and especially biological processes 
have not been addressed in detail. A realistic model scheme ("the probabilistic two- 
stage model") based on probabilistic description of the main processes involved in the 
corresponding radiobiological mechanism will be presented in this paper. 

The probabilistic two-stage model is based, in principle, on distinguishing the single- 
particle and collective effects of the radiobiological mechanism: Individual events of 
energy transfer, subsequent chemical reactions and also formation of DNA damages 
occur practically independently after the impacts of individual ionizing particles, at 
least at dose levels and dose rates that are used commonly in radiotherapy. On the 
other hand, the final response of a given cell, which includes biological processes that 
start from DNA repair and lead finally to cell survival or inactivation, has to be classified 
as a response to the total effect of all traversing particles. 

In the case of protons and other ions, the amount of energy deposited in a cell 
nucleus by a single particle (including all secondary particles) is much higher than 
in the case of irradiating by photon or electron beams. The number of particles that 
contribute to a given dose level is lower, and the severity of damages formed by individual 
particles increases significantly. Even though it is still the collective effect of all the 
damages that leads to cell inactivation, the role of single-particle effects is much more 
important than in the case of photon irradiatio n. Even the impac t of a single particle 
may result in inactivating the cell; compare e.g. ((Miller et al 1119991 ) . Model scheme that 
enables to represent the effects of single particles as well as their synergetic or saturation 
combinations is very helpful in such a case. 

The basic idea of distinguishing the single-particle and collective effects in the 
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case o f proton and light ion irradiation has been formulated by lJudas and Lokaifcek 
(2001), who also suggested that survival curves should be represented by higher-order 
polynomials instead of their simple parabolic description by the LQ model. Further 
analyses of the model framework have demonstrated that the polynomial expansion can 
be used reliably in the low-dose region only The model has been reformulated in order 
to ensure its applicability to survival cu rves measured over a wide range of doses; it 
has been also developed to greater detail (|Kundratl l2004). Applications of the model to 
analyses of experimental data have shown the possibility to describe a large variety of 
experimental cell survival curves. In comparison to phenomenological models, the two- 
stage model has enabled to represent systematically not only the global shape of cell 
survival curves but also their detailed structure. The importance of this feature of the 
model arises from the fact that even small deviations of cell survival curves from their 
global parabolic or linear shape might be largely amplified in the case of fractionated 
irradiation. 



2. Probabilistic two-stage model 

2.1. Basic model framework 

Within the probabilistic two-stage model of biological effects of ionizing particles, the 
following characteristics of the underlying radiobiological mechanism are taken into 
account: (i) the actual number of primary particles depositing energy to the cell nucleus, 
(ii) the amount of energy deposited (including, for the sake of simplicity, the energy 
delivered by all secondary particles), and (iii) the resulting probability of cell inactivation 
which reflects the complexity of the DNA damage produced. 



2.1.1. Number of traversing particles. The numbers of particles traversing nuclei of 
individual cells of the irradiated tissue or sample are of stochastic nature. Their average 
number, k av , is given by particle fluence $ or applied dose D, 

k av = hD . (1) 

Here, h denotes the proportionality constant, which reflects also geometrical 
characteristics of the given cell line and irradiation arrangement. The distribution of 
traversing particle numbers over the cell population can be described then by Poisson 
statistics, 

P k (D) = ¥^eM~hD) . (2) 

Here we have assumed that the impacts of individual particles can be considered as fully 
random, and that the diameters of the most efficient track parts are by far smaller than 
the size of cell nucle us, which is in agreement with the results of track structure studies 
faollev et allll990h . 
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2.1.2. Energy transferred to cell nucleus. The amount of energy transferred to cell 
nucleus (or, more specifically, to chromosomal system and sensitive regions within it) 
together with the spatial distribution of individual energy transfer events play key roles 
with respect to the complexity of the resulting DNA damage. This is in turn crucial 
for the future fate of the given cell. The amount of energy transferred to chromosomal 
system by each primary particle (including all its secondaries) is a stochastic quantity. 
It is influenced by energy and energy-loss spectra of beam particles. It may be, in 
principle, estimated on the basis of theoretical predictions, Monte Carlo calculations 
or microdosimetry measurements. Let us denote by 7r 1 (e) the spectrum of energy e 
deposited to cell chromosomal system in one event (one traversal, i.e. one primary 
particle plus all its secondaries). The normalization condition reads 

lax 

de7n(e) = l, (3) 

where e max stands for the maximum energy that may be transferred in one event. The 
spectrum of energy transferred in k events is given, then, by a corresponding convolution 



7Ti * 7Ti * . . . * 7Ti 

V v ' 

k 

ch Sr, 



/ de' 7n(eV fc _i(e-e') , (4) 
Jo 



which fulfils the condition f £max de irk(e) = 1. The average energy deposited per cell 
nucleus (or chromosomal system), which is obtained by weighting correspondingly the 
average energy deposited by a given number of particles, 



fk Smax 

£ av = Y,Pk / de e ir k (e) , (5) 
k Jo 

is proportional to the applied dose D (or to particle fluence $). 



2.1.3. Cell survival probability. The number of events (traversals of primary particles), 
k, and the total amount of transferred energy, e, stand for the primary characteristics 
describing processes of the physical phase of the radiobiological mechanism. The 
subsequent chemical and biological processes include large cascades of chemical reactions 
that lead to the formation of DNA damage, DNA repair processes, and further 
processes of the biological response of the cell that result finally either in its survival 
or inactivation. As far as cell survival is concerned, these complex processes can be 
taken into account in a simplified way by representing the final endpoints directly, i.e. 
by considering only the final probabilities of cell inactivation. In the following we shall 
denote by (e) the average cell inactivation probability after the impact of k (primary) 
particles that have deposited energy e. Note that a similar approach, based e.g. on DSB- 
induction probabilities, might be used if other effects, e.g. the yields of DNA damage, 
were the endpoints of interest. Probability of cell survival (i.e., cell survival curve) is, 
then, given by 

fk Smax / -\ 

s(D) = l-J2Pk(D) de n k (e)p$\e), (6) 

k Jo 
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where the distribution of the actual number of particle traversals, k, and also the 
spectra of transferred energy, 7Tfc(e), have been taken into account. Here the dose 
dependence is included in Pk(D), see Eq. (j2J). Note that the average number of 
traversing particles per unit dose, h, as well as transferred energy spectra, vr fc (e), and 
also inactivation probabilities (e) vary significantly for different particles, their LET 
values and/or diverse cell lines. This is in agreement with the experimental evidence of 
the dependence of cell survival on all these factors. Given the inactivation probabilities 
Pk\e) together with the distribution of particle numbers Pk and the spectra of 
transferred energy 7Tfe(e) under the given experimental or clinical irradiation conditions, 
one would be able to calculate the c ell survival probability according to Eq. © 



Using corresponding model a pproaches ((Webb and Nahumlll993llNiemierko et al 



1992 



Niemierko and Goiteinl I1993T ). macroscopic endpoints such as the tumor control and 



normal tissue complication probabilities (TCP and NTCP) can be estimated, then. In 
evaluating the cell survival probability according to Eq. (JBJ), a crucial role is played by 
the inactivation probabilities p®. These probabilities reflect the biological effectiveness 
of different radiation modalities. They may be interpreted directly in terms of physical 
and biological processes. They represent, therefore, a solid basis for establishing more 
detailed microscopic models of radiobiological mechanism. Methods of evaluating these 
probabilities on the basis of analyzing corresponding experimental cell-survival data will 
be derived and discussed in the next sections. 



2.2. Monoenergetic ions 

Radiobiological experiments aimed at establishing the biological effectiveness of different 
radiation kinds are performed usually by irradiating cell monolayers by monoenergetic 
beams only. The beam energy is chosen so as to correspond to certain regions in the 
Bragg peak of the given particles (so-called track-segment experiments). In such a 
case the spectrum of transferred energy is rather narrow, at least as compared to non- 
monoenergetic beams. The (average) energy transferred to cell nuclei is proportional to 
the LET value A of the particles. Moreover, in this case the average number of (primary) 
particles traversing cell nuclei, i.e. the parameter h (cf. Eqs. ((H) and (J2J)), can be related 
to the effective cross-section, a, of the cell nucleus (or chromosomal system) and LET 
value A by 

Approximating the density of traversed medium by that of water, the conversion 
constant C = 6.24 keV Gy _1 /iin~ 3 for the usual choice of units (A[keV//im], _D[Gy], 
cr[/im 2 ]). Note that a denotes the geometrical effective cross-section of the chromosomal 
system (or of sensitive region within cell nucleus) to be traversed by the given particles. 
The consequent processes of DNA damage and cellular response are not involved in 
a, but treated separately. Note especially that o is not identical to cell inactivation 
cross-section, which is used by some authors to describe (the linear component of) the 



Probabilistic two-stage model of cell inactivation by ionizing particles 



6 



proba bility, per 1 particle, of the cell to be inactivated; compare e.g. (|Scholz and Kraft 



1994). 



Cell survival probability can be, then, expressed as 

siP) = p kP, A) (l -pf>(A)) = 1 - E^(A A) pj»(A) , (8) 



where (average) inactivation probabilities p& (A) are now expressed as functions of LET 
value A instead of energy e; they are given by pi (A) = J^ £max de 7Ck{e) p k 9(s). 



2.3. Detailed model of inactivation probabilities 

In the case of particles characterized by higher LET values, i.e. protons, neutrons and 
ions, even the impact of a single primary particle (including the effects of all its secondary 
particles) may produce a severe damage to chromosomal DNA, which might finally lead 
to cell inactivation. In the following, we shall refer to these damages as to "single- 
particle induced" damages. Let us denote the probability of their induction, per one 
traversing particle, by a. If k particles have traversed the chromosomal system, the 
probability of not being damaged in this way is given by 

q£ = (1 - a) k . (9) 

A vast majority of DNA damages formed by traversing particles are not of this 
severity; they can be repaired successfully and do not lead to cell inactivation. However, 
DNA damages caused by different particles might combine to form a kind of damage 
that is difficult to be repaired. Synergetic combinations of such damages might include 
combinations of single-strand breaks (SSBs) into a double-strand break (DSB) or DSBs 
into a more complex damage. They need not be limited to one chromosome only; e.g. the 
impact of a particle pair damaging corresponding segments of chromosomes of a given 
homologous pair. To account for these combined effects, let us denote by b the average 
probability that such a severe damage, which alone does not result in inactivating the 
cell, has been formed by a single particle. If this is the only damage to the cell, we shall 
assume that it can be repaired quite easily and the cell will survive, possibly with some 
delay in its proliferation cycle. However, combination of two or more such damages 
might lead to cell inactivation if not repaired correctly. Probability that no such "pair" 
or "combined" damage^ has been formed after the impact of k particles is then, to the 
first approximation, given by 

q% = (i _ &2)*(fc-D/2 . (10) 

Here we have assumed, for the sake of simplicity, that the synergetic effects occur 
predominantly from particle pairs only; the exponent in the last formula stands just 
for the number of particle pairs. 

I Note that also the probability of inducing "pair damage" to DNA refers to the effect of a single 
particle. The term "pair" or "combined" damage has been chosen just to highlight the lower severity of 
such damage, i.e. the necessity of being combined with another damage of a similar character in order 
to lead to cell inactivation. 
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The final biological effect, measured e.g. by means of cell survival, does not depend 
only on the initial yield of DNA damages. Cellular repair processes have been shown to 
contribute significantly to this point, too. In the following we shall assume that severe 
DNA damages potentially leading to cell inactivation, either caused by single particles 
("single-particle induced" damages) or combined from multi-particle effects ("combined" 
or "pair" damages), may be repaired by the cell. We shall denote the probabilities of 
repairing them successfully by r A and r B . After the cellular repair processes took place, 
the probability that a residual severe damage is present equals, then, 

Pt = (1 - §?)(! - r£) = (l - (1 - a) k ) (l - r£) (11) 

in the case of single-particle induced damages, and 

p B k = (1 - §?)(! - rf ) = (l - (1 - b^~^) (l - rf ) (12) 

in the case of combined (or "pair") damages, respectively. Here we have explicitly 
indicated that the rate and fidelity of repair processes in the two cases need not be 
equal. The damage complexity increases with the number of traversing particles k and 
with the LET value A. Therefore, the probability of successful repair of DNA damages 
should decrease with these factors. The repair probability might, however, in some cases 
exhibit more complex behavior, corresponding to the onset of cellular repair processes. 

The probabilities that the mentioned damages have either not been formed at all 
or have been repaired faithfully are just the complements of the last terms, 

^ = l-^ = l-(l-(l" <»)*) (l " r£) (13) 

and 

q£ = 1 - p* = 1 - (l - (1 - br ik - 1)/2 ) (l - rf ) , (14) 

respectively. If the cell should survive, it must not be inactivated by either of the 
mentioned mechanisms; i.e., neither single-particle induced nor pair damages must be 
present after the repair processes took place. The cell survival probability is, therefore, 
given by 

% = qt q% , (15) 
and the inactivation probability 

pf = 1-Qk- (16) 

The cell survival probability s(D) after applied dose D is given then by Eq. (JHJ), 
where the distribution of the numbers of particles traversing individual cell nuclei is 
given by Poisson distribution, Eq. (J2J). In the case of Bragg peak or other non- 
monoenergetic irradiation, the spectra of transferred energy have to be taken into 
account, too, according to Eq. (ffi|) . Similarly, in the case of mixed beam irradiation 
one has to convolve these results with the spectrum of traversing particles (i.e. with the 
abundances of individual ion species in the beam), again in the way analogous to that 
indicated in Eq. (JHJ). 
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The probabilities of "single-particle" and "pair" damage formation, a and 6, as well 
as of successful repair, and r^, depend significantly on the radiation quality (ion 
kind and LET value) and, of course, on the biological characteristics of given cell line 
(cellular radiosensitivity, cell cycle phase etc.). The values of these parameters may be 
established by analyzing corresponding experimental data. Systematic sets of data, in 
which cell survival has been assessed for a given ion at different energies (different LET 
values), might be used to evaluate even the functional dependence of these probabilities 
on particles' LET values. Corresponding analyses have been performed by us for various 
ions; the results will be shown in the following sections. In determining the values 
of the damage induction and repair probabilities, more detailed microscopic models, 
representing the underlying chemical and biological processes, might be helpful, too. In 
fact, the given model scheme represents a framework enabling one to relate the outputs 
of those microscopic models to cell survival characteristics, which might be assessed 
experimentally in a direct way. 

The probability of successful repair, r, should reflect the complexity of the total 
damage formed by all the traversing particles. The damage complexity might be, to 
the first approximation, estimated from the total amount of energy deposited to cell 
nucleus (or chromosomal system). For the sake of simplicity, we have therefore limited 
ourselves to consider the repair probability, r(k, A), to be a function of the transferred 
energy only, e ~ k.X: 

r(k, A) — > r(k.X) . 



2.4- Basic model features 

Let us discuss the basic features of the given model framework, especially the role 
of individual model parameters in determining the shapes of cell survival curves. 
Obviously, there are no difficulties in describing simple linear or parabolic cell survival 
curves by the given model scheme. Linear curves (Fig. P) can be represented by taking 
into account the terms including p£ only, i.e. by assuming that all the lesions formed 
by the traversing particles are lethal, do not need to be combined with damages caused 
by other particles (i.e. 6 = 0), and cannot be repaired by the cell (r = 0). This 
is in agreement with experimental evidence: Linear survival curves have been found 
especially in very sensitive cel l lines and in cell lines that have been proved to be 
deficient in repair (j,Ioinerlll997h . Linear response has been observed also in irradiating 
by ve ry high-LET parti cles, where the severity of DNA damages is assumed to be very 



high (jKraft et a/lll997f ). Similarly, survival curves characterized by a global parabolic 
shape can be represented by taking into account the terms including pf ; compare Fig. El 
Even cell lines exhibiting DNA repair of a certain level may be characterized by roughly 
parabolic survival curves; compare Fig. El 

Cell survival curves may, however, exhibit more complex behavior. The given model 
scheme can be used to represent survival curves showing, e.g., a downward bending in 
the region of lower doses followed by an upward bending at higher doses, compare 
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Figure 1. Examples of linear cell survival curves, obtained by the probabilistic two- 
stage model by putting h — 2; a = 0.2 (solid line) or 0.5 (dashed line); 6 = 0; 
r A = = r B . Note that the scaling of the cell survival with applied dose is somehow 
arbitrary in these examples; it corresponds to the choice of h value. 



□ [Gy] 
2 1 " 3 




0.01 J 



Figure 2. Parabolic cell survival curves corresponding to h = 2; a = 0.1; b = 0.3 
(solid line) or 0.8 (dashed line); repair not taken into account, r A = = r B . 



Fig. |3J The corresponding repair probability, to which such complex behavior has to be 
attributed, is shown, too. 

The model enables one to describe also survival curves exhibiting the low-dose 
hypersensitivity phen omenon, which has been observed for many cell lines irradiate d 
by different particles (j,Ioiner et al 1120011 IChalmers et al 112004 ISchettino et al 1120011 ). 
In Fig. |S] two examples of calculated survival curves with hypersensitivity regions are 
demonstrated. The corresponding repair characteristics, showing an onset (triggering) 
effect, are displayed, too. Note that if the repair profile exhibits significant variations in 
the rather narrow range of a few particles traversing the nucleus, the corresponding 
survival curve may be characterized even by a temporary increase in the medium- 
dose regfon_(FigJ^]b Such behavior has been found in several experiments, compare 
e.g. ([Chalmers et al Ibooi. 

The examples shown in this section demonstrate clearly that in determining the 
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Figure 3. An example of survival curve characterized by a global parabolic shape 
even if repair taken into account. Damage induction probabilities a = 0.1; b = 0.8; 
number of particles traversing the nucleus per unit dose h = 2. Repair of combined 
damages as shown in the right-hand panel; parameterized by auxiliary function 
r{k) = cxp(— (rife)*" 2 ), where n = 0.2, r 2 = 2. Repair of single-particle induced 
damages not taken into account, r A = 0. 
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Figure 4. Example of more complex behavior of cell survival curves. Solid line 
obtained by putting h = 1, a = 0.1, b — 0.8. Repair probability, decreasing with 
the number of particles traversing the nucleus, as depicted in the right-hand panel; 
parameterized by auxiliary function r{k) — exp(— (rifc) r2 ), where r\ — 0.5, r 2 = 0.4. 
Survival curve shown by dashed line corresponding to the same parameter values but 
without repair, r B — 0; in both the cases being put r A — 0. 



shape of cell survival curves the repair characteristics may be even more important than 
the damage induction itself. 
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Figure 5. Examples of survival curves exhibiting low-dose hypersensitivity 
phenomenon. Here h = 4, a = 0.4, 6 = (solid line) and h = 6, a = 0.5, b = 0.2 
(dashed line), respectively. The repair probability r A (right panel) is characterized by 
onset behavior, with maximal repair rate in the region of a few traversing particles. 
Parameterized by auxiliary function r(k) = exp(— (rifc) r2 )[l — cxp(— (r 3 fc) r4 )], where 
r = 2, n = 0.2, r 2 = 2, r 3 = 0.3, r± = 3 (solid line) or r = 1, r\ = 0.15, r 2 = 1.5, 
r 3 = 0.5, r4 = 5 (dashed line); r B = in both the cases. 




Figure 6. Example of model representation of survival curve exhibiting highly 
pronounced hypersensitive behavior, with a temporary increase following the low-dose 
hypersensitivity region. Probability of single-particle induced DNA damage a = 0.5, 
average number of traversing particles per unit dose h — 6. Combined damages not 
taken into account, b — 0, r B — 0. Auxiliary parameterization of the repair probability 
(right panel) r(k) = cxp(— (rifc) r2 )[l — cxp(— (r 3 fc) r4 )], where r = 1, r\ = 0.1, r 2 = 2, 
r 3 = 0.5, r 4 = 3. 

3. Application to experimental data 

The proposed model has been applied to analysis of experimental cell survival data for 
Chinese hamster V79 cells after irradiation by low-energy protons, which was assessed 
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experimentally bv iBelli et al I ([19981 ). Mono-energetic proton beams were used in their 
study, with the energies in the range of 0.57 - 5 MeV, corresponding to LET values 
7.7 - 37.8 keV//im. 

The following parameterization has been used for representing the damage induction 
and repair probabilities: 

a(X) = (aiX + a 2 A 2 ) 1 - exp ( - (a 3 A) a4 ) , 

b(\) - 1 ~ ex PHM) b2 ) {m 

6(A) " l + 6 3 exp(-(MN ' (17) 
r*(Jb,A)- ^^-W) 



1 + r 3 exp( — (r 4 fcA) r5 ) 

Flexible monotonous test functions were used in 6(A) and r B (k,X), whereas a(A) has 
been allowed to possess a non-monotonous behavior. 

Model calculations of the survival curves (full lines) together with experimental data 
(points) are shown in Fig. [7| For comparison, fits of the data according to the linear- 
quadratic (LQ) model are displayed, too. Probabilities of forming "single-particle" and 
"pair" DNA damages discussed in the preceding are shown in Fig. |HJ Single-particle 
induced damages, a, have turned out to play a minor role only. Consequently, it has 
been possible to neglect their repair, i.e. to put r A (k, A) = 0, which has simplified the 
fitting procedure. The derived probability of successful repair of the "pair damages" as 
a function of traversing particle numbers k and LET values A is plotted in Fig. |H1 

The optimization procedure used in the present analysis has been based on 
applying standard minimization methods; mainly t he SIMPLEX and MIGRAD methods 
implemented in the MINUIT minimization tool ((Jame s 1994) have been applied to 
dedicated computer codes written in the programming lang uage FORTRAN. The 
goodness of fit has been assessed by calculating the \ 2 values (jHagiwara et al 1 12002) 
for individual survival curves. They have been calculated as a sum of the deviations 
of (logarithms of) model calculations, s mo d e i, from the data, s exp , weighted by the 
experimental errors, A,, of individual data points (again in the logarithmic scale), 

2 \ ^ (111 Smodei hi S eX p) 

i i 

The values of \ 2 for individual survival curves represented by the probabilistic two-stage 
(P2S) model and by the linear-quadratic (LQ) model are listed in Table The total 
X 2 values were 24.0 for the P2S model and 45.0 for the LQ model, respectively. The 
corresponding values of auxiliary parameters involved in Eqs. ()17j) are listed in Table 
These values correspond to the effective cross-section of V79 nucleus (or chromosomal 
system) a = 12.8 /im 2 ; the goodness of fit of the model to the data has varied only 
slightly around this value. 
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Figure 7. Survival curves fo r V79 cell s irrad iated by low-energy protons. 
Experimental data established by Belli et al \ (^998), interpreted on the basis of a 
systematic analysis using the detailed inactivation scheme of the probabilistic two- 
stage model (left panel). The corresponding probabilities of DNA damage formation 
and repair are shown in Figure [SJ Fits by the LQ model shown, too (right panel). 



Table 1. Experimental survival curves and precision of representing them by the 
probabilistic two-stage (P2S) and the linear-quadratic (LQ) models. 



Energy [MeV] 


5.01 


3.20 


1.41 


0.76 


0.64 


0.57 


LET [keV/^m] 


7.7 


11.0 


20.0 


30.5 


34.6 


37.8 


XP2S 


2.78 


3.84 


0.70 


9.53 


4.44 


2.73 


Xlq 


2.98 


4.36 


1.08 


6.06 


14.4 


15.8 
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Figure 8. Left panel: The probabilities of single-particle induced, a (solid line), 
and pair DNA damages, b, as functions of proton LET values. Right panel: The 
probability of successful repair, decreasing with both LET and number of traversing 
protons, reflecting the increasing complexity of DNA damage. 



Table 2. Values of auxiliary parameters involved in model representation of damage 
induction and repair probabilities, Eq. Q17|). 



Single-particle induced damage 


Pair damage 


Repair probability 


o(A) 


6(A) 


r B (k,\) 


a x = 0.0022 /Ltm/keV 


6i = 0.12 fim/keV 


ri = 0.024 /zm/keV 


a 2 = 0.013 (^m/keV) 2 


6 2 = 5.0 


r 2 = 5.0 


a 3 = 0.026 /zm/keV 


6 3 = 24.1 


r 3 = 0.56 


<24 = 5.0 


64 = 0.061 fim/keV 


r 4 = 0.0021 /im/keV 




65 = 1.76 


r 5 = 5.0 



4. Discusion 



Basic characteristics of DNA damage formation and of cellular repair processes have 
been derived in the studied case. The corresponding probabilities have been determined 
by analyzing experimental cell survival data only; no other inputs or assumptions have 
been used. Nevertheless, the results are co nsistent with expe rimental studies concerning 
the yields of DNA damage; compare e.g. ()Prise et al Ill998f ). 

The results have shown that for protons in the studied range of rather high LET, the 
pair damages dominate over the single-particle ones. In the local-effect model (LEM), for 
example, combined effects are not taken into account at all, as all damages are assumed 
to be of lethal nature. This seems to be the reason why discrepancies of the order 
of 25% have been found between the given experi mental data and t he corresponding 
calcu lations based on the LEM scheme; compare (jBelli et al I Il998l IScholz and Kraft 
19951 ) . The differences in inactivation mechanism between protons and heavier ions, 
identified on the basi s of the given detailed model scheme, will be discussed in detail in 
a forthcoming paper (jKundrat et al 1 12004). 

The formation of lethal damages by individual particles (single-particle induced 
damages) gets saturated for protons at LET values around 30 keV//im. This means 
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that for protons the RBE effects related to overkill occur already on the level of 
damage induction probability, which has a maximum at values much lower than unity 
(around 20% in the present analysis). The probability of inducing combined damages, 
on the other hand, possesses a monotonous increase over the whole studied LET 
range. The repair probability exhibits decreasing behavior with increasing LET and/or 
particle numbers, corresponding to increasing complexity of the total damage induced 
in chromosomal DNA. 

The effective cross-section of cell nuclei has been established in the present analysis 
on a phenomenological basis only. It has been limited from below by the given data, 
as by assuming a too low value of the cross section the cell survival probability would 
be too high as compared to the given data. Upper limit on the effective cross section 
has been represented by the average geometrical cross-section o f cell nuclei of the given 
cell line a geom = 134keV//im reported in the given experiment ( Belli et al 199$ ). The 
effective cross-section (i.e., the cross-section of the sensitive region within the nucleus 
with respect to proton tracks) has corresponded to approximately 10% of the nucleus. 

In the future, attempts will be made to relate the damage induction and repair 
probabilities to the results of track structure studies and other microscopic models of 
physical and biological processes. These models might be helpful in determining the 
effective cross section of the cell nuclei with respect to a given ion (and its energy), too. 
The present analyses rely on the fact that the most efficient parts of proton tracks are 
much narrower than the dimensions of cell nuclei, which has enabled to apply Poisson 
statistics. 

As to the precision of representing experimental cell-survival data, the present 
calculations and those based on the LQ model yield similar results for lower-LET 
protons, where the experimental survival curves possess quite simple shapes. They 
differ significantly in the highest LET values (compare Fig. |7J) corresponding to proton 
track ends, where detailed structure of the data can be represented only by the detailed 
scheme of the probabilistic two-stage model. 

The required level of accuracy influences the number of auxiliary parameters that 
have to be involved in the parameterization of DNA damage induction and repair 
probabilities. If e.g. in treatment planning applications a less precise representation 
of survival curves is sufficient, less detailed parameterization e.g. of the form of 
f(x) = exp(— (ax) 13 ) might be used for 6(A) and r B (k, A) instead of those of Eqs. (fT7|) . 
This difference in parameterization is of technical nature only, the basic scheme of the 
model remaining the same. 



5. Conclusion 



The probabilistic two-stage model provides a realistic description of the radiobiological 
mechanism. It includes description of DNA damages of different severity caused by 
individual traversing particles, as well as of their synergetic and saturation combinations. 
Basic characteristic of repair processes are incorporated, too. The model might serve as 
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a basis for more detailed microscopic modeling of radiobiological effects. 

As demonstrated explicitly for the case of low-energy protons, the probabilistic 
two-stage model enables to represent experimental cell survival curves with much 
better precision than the LQ model, while the total number of parameters involved 
is about the same. The ability of representing cell survival precisely is a prerequisite 
if different fractionation schemes of hadron radiotherapy are to be evaluated, because 
in fractionated irradiation the local variations (i.e. local deviations of the cell survival 
ratio from the global shape of survival curve) may get largely amplified in the cumulative 
effect of the given scheme. The model enables to represent simple linear and parabolic, 
but also rather complex survival curves, including e.g. the low-dose hypersensitivity 
phenomenon, where a kind of triggering behavior of repair probability has to be taken 
into account. Unconventional fractionation schemes including hyperfractionation might 
be assessed on this model basis, too. 
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